arXiv:1502.05944v2 [cond-mat.stat-mech] 30 Apr 2015 


Thermodynamic Uncertainty Relation for Biomolecular Processes 


Andre C. Barato and Udo Seifert 

II. Institut fur Theoretische Physik, Universitdt Stuttgart, 70550 Stuttgart, Germany 

Biomolecular systems like molecular motors or pumps, transcription and translation machinery, 
and other enzymatic reactions can be described as Markov processes on a suitable network. We 
show quite generally that in a steady state the dispersion of observables like the number of con¬ 
sumed/produced molecules or the number of steps of a motor is constrained by the thermodynamic 
cost of generating it. An uncertainty e requires at least a cost of 2fcsr/e^ independent of the time 
required to generate the output. 
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Biomolecular processes are generally out of equilibrium 
and dissipative, with the associated free energy consump¬ 
tion coming most commonly from adenosine triphosphate 
(ATP) hydrolysis. The role of energy dissipation in a va¬ 
riety of processes related to biological information pro¬ 
cessing has received much attention recently , to 

give just one class of examples for which one tries to un¬ 
cover fundamental limits involving energy dissipation in 
biomolecular systems. 

Chemical reactions catalyzed by enzymes are of cen¬ 
tral importance for many cellular processes. Prominent 
examples are molecular motors |15l - ll9l |. which convert 
chemical free energy from ATP into mechanical work. 
In this case an observable of interest is the number of 
steps the motor made. Another commonly analyzed 
output in enzymatic kinetics is the number of product 
molecules generated by an enzymatic reaction, for which 
the Michaelis-Menten scheme provides a paradigmatic 
case 0. 

Quite generally, chemical reactions are well described 
by stochastic processes. An observable, like the rate of 
consumed substrate molecules or the number of steps of 
a motor on a track, is a random variable subjected to 
thermal fluctuations. Single molecule experiments na¬ 
il^ provide detailed quantitative data on such random 
quantities. Obtaining information about the underlying 
chemical reaction scheme through the measurement of 
fluctuations constitutes a field called statistical kinetics 
[25l - [^ . A central result in this field is the fact that 
the Fano factor quantifying fluctuations provides a lower 
bound on the number of states involved in an enzymatic 
cycle [TsI. [^. 

For a non-zero mean output, the chemical potential dif¬ 
ference (or affinity) driving an enzymatic reaction must 
also be non-zero, leading to a free energy cost. Is there a 
fundamental relation between the relative uncertainty as¬ 
sociated with the observable quantifying the output and 
the free energy cost of sustaining the biomolecular pro¬ 
cess generating it? 

In this Letter, we show that such a general bound does 
indeed exist. Specifically, for any process running for a 
time t, we show that the product of the average dissi¬ 
pated heat and the square of the relative uncertainty of 


a generic observable is independent of t and bounded by 
2kBT. This uncertainty relation is valid for general net¬ 
works and can be proved within linear response theory. 
Beyond linear response theory, we show it analytically 
for unicyclic networks and verify it numerically for sev¬ 
eral multicyclic networks. As an illustration of a specific 
consequence of our results, we obtain a new bound on 
the Fano factor for unicyclic networks which depends on 
the affinity driving the process. 

The observables we consider here arise from counting 
the number of transitions between states, which is differ¬ 
ent from a random variable counting the fraction of time 
spent in a state. For this latter quantity, the relative un¬ 
certainty is finite even in equilibrium, as is the case of 
a cell estimating the concentration of an external envi¬ 
ronment by counting the fraction of time receptors are 
bound [1^, [1^ . The role of dissipation for this problem 
has been recently studied in [l^ . 

Our main result can be introduced using the arguably 
simplest example of a nonequilibrium chemical reaction 
catalyzed by an enzyme, which is a biased random walk 
where a single step is interpreted as the completion of 
an enzymatic cycle. Steps to the right happen with a 
rate k'^, those to the left with a rate k~. After a time t, 
on average, {X) = (fc+ — k~)t steps have occurred, with 
the number of steps X corresponding to the observable 
of interest. Specific realizations of this rando m p rocess 
show a variance ((A — (A))^) = (fc+ -|- k~)t [sj. The 
squared relative uncertainty of the observable is 

= ((A2)-(A)2)/(A)2 = {k+ + k-)/[{k+-k-)H]. (1) 

Assuming an external environment of fixed temperature 
T, the transition rates are given by the local detailed 
balance relation 

k'^/k~ = eyi.-p{A/kBT), (2) 

where A is the affinity driving the process and kB Boltz¬ 
mann’s constant. The thermodynamic cost of generating 
this output is given by the entropy production rate, which 
reads 

a = {k+- k-)A/T, (3) 

leading to a total dissipation after time t of Tat. In 
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equilibrium, i.e., for k~^ = k , there is no dissipation and 
the uncertainty e diverges. 

The trade-off between precision and dissipation is cap¬ 
tured by the crucial product Q of total dissipation and 
the square of the relative uncertainty, 

Q = Tate^ = A coth[A/(2kBT)] > 2/cbT, (4) 

where we used Eqs. o, ©, and ©. This thermody¬ 
namic uncertainty relation shows that a more precise out¬ 
put requires a higher thermodynamic cost independent of 
the time used to produce the output. Reaching an un¬ 
certainty of, e.g., one percent requires at least 20000 A:_bT 
of free energy. Since Q is an increasing function of the 
affinity A, the minimal cost for a given uncertainty is 
achieved close to equilibrium, i.e., for ^ 0. In this 

limit, however, the time t required for producing a sub¬ 
stantial output {X) diverges. 

In the following, we show that this uncertainty relation, 
namely, that the dissipation of a process that leads to an 
uncertainty e must be at least 2kBT/e^, is quite general, 
holding true for arbitrary networks of states. First, we 
prove Q > 2k bT for any network within linear response 
theory. For a unicyclic network we show analytically that 
the bound also holds true beyond linear response. For 
multicyclic networks beyond the linear response regime, 
we provide numerical evidence for this bound. From now 
on, to keep notation slim, we set /cs = T = 1, which 
renders entropy and energy dimensionless. 

We consider a general Markov process with transition 
rate from state i to j denoted by kij. Thermodynamic 
consistency requires that if kij ^ 0 then kji ^ 0. Fur¬ 
thermore, we assume a finite number of states N and 
denote the stationary probability of state * by P^. 

The observable of interest represents some phys¬ 
ical quantity that changes if certain transitions in the 
network of states occur. Specifically, the generalized dis¬ 
tance determines how much the variable X^ changes 
if the transition i to j happens. As an example, if X^ 
counts the number of consumed ATP molecules, and if 
state i represents a free enzyme and j an enzyme with 
ATP bound to it, then dg = I and = —I. This 
generalized distance is always antisymmetric in i and j. 

The affinity associated with the variable X^ is denoted 
Aa- For example, if X^ is the number of consumed sub¬ 
strate molecules (like ATP) in a chemical reaction, then 
Aa is the chemical potential difference driving this reac¬ 
tion. The transition rates fulfill the generalized detailed 
balance relation [s^ 

\n{hj/kji) = ^ 4jAp +Ei- Ej , (5) 

P 

where this sum is over all affinities Ap , including the case 
(3 = a, and Ei is the equilibrium free energy of state i. 

In the stationary state, the velocity (or probability cur¬ 
rent) and diffusion constant associated with Xa are de¬ 


fined as 

Ja = {Xa)/t (6) 

and 

Da ^ [{XD - {Xaf]/{2t), (7) 

respectively. The squared relative uncertainty then reads 

ei ^ [{XD - {XafyiXa)^ = 2Da/{Jlt). (8) 

While the probability current has a simple form in terms 
of the stationary probability distribution, namely, Ja = 
yZij dfj{Pikij — Pjkji), a general formula for the diffusion 
constant is more involved and will be discussed below. 
The entropy production rate is 

P 

where /3 runs over all affinities. For example, for a molec¬ 
ular motor this sum has two terms: one affinity is the 
chemical potential difference driving the motor with the 
rate of ATP consumption as the associated current, the 
other affinity is the mechanical force and the respective 
current is the velocity of the motor. The (dimensionless) 
product (j4]) for a general network is then defined as 

= atel = 2aDalJl. ( 10 ) 

Within linear response theory , which is valid close 
to equilibrium where the affinities Ap are small, a current 
can be expressed by the affinities as 

dp ~ ^ ^ LfijAj, ( 11 ) 
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where the Onsager coefficients are defined as 

Ljp = = Lpj. ( 12 ) 

From Fqs. (I9|) and m the entropy production within 
linear response reads a = ^ Lp-yA^Ap. Moreover, 

the diffusion constant is given by the Finstein relation 
Da = Laa [i^ , which from Eq. ([8]) leads to being 
linear in the affinities. Hence, Eq. m becomes 

G/SjA^Aj 

'lip-y LapLa-fApA^ V {JaY 

(13) 

where Gp-y = {LaaLpy — LapLay)- Using the fact that 
the Onsager matrix L is positive semi-definite it is possi¬ 
ble to prove that G is also a positive semi-definite matrix 
[s^ . Hence, we have established Qa > 2 within linear re¬ 
sponse theory. Note that equality is reached in the case 
of only one non-zero affinity, i.e., Ap = 0 for /3 a. 

In the calculations that follow we use elegant expres¬ 
sions obtained by Koza [H, for velocity and diffusion 
coefficient, which are valid for a general network of states. 
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For these expressions we need a modified generator asso¬ 
ciated with Xn^ which is a A^-dimensional square matrix 
with elements 




%exp(zdf) 




j "-y 


iii= j 


(14) 


A set of coefficients Cn{z) is defined through the charac¬ 
teristic polynomial of this matrix as 


N 


det {yl - £“(z)) = <^n(z)y", 


(15) 


n—0 


where / represents the identity matrix. Using these co¬ 
efficients, which are functions of the transition rates, the 
current and diffusion coefficient can be written as 13511 




(16) 


and 


Da = -(U" + 2C( J„ + 2C2 J^)/(2Ci), (17) 


where Cn = Cn(0) and the primes denote derivatives 
with respect to z taken at z = 0. A full derivation for 
these expressions is given in 341. 


We first consider an arbitrary uni-cyclic model with X 
states [15 . The transition rate from state n to state n-l-1 
(n — l) is denoted (fc“), where n = 0,1,..., — 1. 

The output X counts the number of completed cycles. It 
is sufficient to count the number of transitions through 
one of the links in the cycle, which we choose to be the 
link between states 0 and 1. The generalized distance 
associated with X is then c?oi = —dio = 1 and dij = 0 
for ij ^ 01. The cycle affinity is 


A = in(r+/r_), (18) 

where r+ = and T" = ^i- example 

of such unicyclic machine with iV = 3 is an enzyme E 
that consumes ATP according to the scheme 

E+ATP ^ ET ^ ED+P^ ^ E+ADP+Pi, (19) 

where ADP stands for adenosine diphosphate, Pi for 
phosphate, ET {ED) represents the enzyme with an 
ATP (ADP) bound to it. In this case the variable X is 
the number of consumed ATP molecules and the affin¬ 
ity is given by the chemical potential difference A = 
y-ATP — y'ADP — yp- 

We can show that for a given number of states N and 
affinity A the product Q reaches its minimal value for 
uniform rates, i.e., kf = (r+)i/^ and fc~ = (r“)i/^ 
independent of f, leading to the bound [sj 

Q > (A/N) coth[A/{2N)] > 2, (20) 

which is in agreement with Eq. (|3]) that corresponds to 
TV = 1. The bound (PIB gives the minimal dissipation re¬ 
quired to realize an uncertainty e for given affinity A and 


number of states N. This bound is an increasing func¬ 
tion of A, hence Q is minimal for A —t 0 where Q —5> 2. 
A related quantity, defined as the ratio of a “baromet¬ 
ric” force and an “Einstein” force, has been considered 
in [15, where a bound similar to Q > 2 has been shown 
to hold for the case N = 2 within a calculation keeping 
terms up to second order in the affinity A. 

We now turn to a specific example showing how this 
new constraint involving fluctuations and energetic cost 
can be turned into a diagnostic tool to unveil a struc¬ 
tural property of an enzymatic cycle. A quantity closely 
related to the relative uncertainty is the Eano factor 

F = [(A^) - (A)2]/(A) = 2D/J, (21) 

which gives a measure of the dispersion associated with 
X, where X counts the output of an enzymatic cycle. 
Eor unicyclic networks, this Fano factor is known to be 
bounded from below by 1/N [H, [l^l- Measurements of 
the Fano factor can then be used to obtain a bound on 
the number of states of an underlying enzymatic cycle 
[ 15 . Our new bound (pH)) implies 

F = QjA > (1/A) coth[A/(2A)]. (22) 

For a diverging affinity, which is the case in chemical 
reaction schemes where at least one backward transition 
rate is assumed to be zero, this bound becomes the known 
one F > 1/N. For experiments where substrate and 
product concentrations are kept fixed and, consequently, 
the value of the affinity is known, as for example in [40l| . 
our stronger bound in Eq. (j22l) constrains even further 
the number of states in such an enzymatic cycle. 

Let us turn again to multicyclic networks. Within the 
linear response regime, we have established above (af¬ 
ter Eq. (HH)) that Qa reaches the bound 2 for the case 
where all affinities but Aa are zero. Eor unicyclic net¬ 
works, which is the paradigmatic case for a system with 
one affinity, we have proved that the bound holds arbi¬ 
trarily far from equilibrium, being reached only in the 
linear response regime. In order to provide full evidence 
that our main result Qa > 2 is indeed universal, we now 
analyze multicyclic networks beyond linear response. In 
this case, we have to take specific systems. As a first 
example we consider a model with an enzyme E that can 
consume two different substrates Si and S 2 and gener¬ 
ates product P, see Eig. [TJ Two enzymatic cycles of this 
model are 

E + Si ^ ESi ^ EP ^ E + P 

k21 ^42 ki4 

E + S2 ^ ES2 ^ EP ^ E + P, (23) 

^31 /C 43 ki4 

where the enzyme states are identified as U = 1 , ESi = 
2, ES 2 = 3 and EP = 4. The affinity of the cy¬ 
cle involving substrate Si {S 2 ) is given by the chemi¬ 
cal potential difference Ai = yi — yp {A 2 = y 2 — yp)- 
The relations between these affinities and the transi¬ 
tion rates are Ai = ^[^ 12 ^ 24 ^ 41 /(^ 21 ^ 42 ^ 14 )] and A 2 = 
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FIG. 1. Multicyclic network of states for an enzyme E that 
can bind substrates Si and S 2 , and produces P. The solid 
(dashed) cycle represents the first (second) cycle in Eq. (1231) . 


ln[fci 3 fc 34 fc 4 i/(A; 3 ifc 43 A:i 4 )]. There is also a third cycle 
1—i'2—>-4—>-3—>-1, in which an is consumed 
and an S 2 produced. Its affinity is not independent but 
rather given by Ai — A 2 - 

The average rate at which 5'i molecules are consumed 
is given by the current Ji = Tifci 2 — ^ 2 ^ 21 , while the rate 
of S 2 consumption is J 2 = Pikis — P^k^i- The entropy 
production can be conveniently written as a sum over 
terms which are the product of a cycle affinity and a 
current (^ . For the present model it becomes 

o = AiJi + A2J2- ( 24 ) 

As an observable of interest we choose the number of 
consumed molecules, for which we can calculate the 
associated product Qi with formulas (nnD, (HU), dni), and 
(l24l) . The resulting function of the ten transition rates 
is too cumbersome to show. Both by minimizing this 
function numerically and by evaluating it for randomly 
chosen transition rates we find that the uncertainty rela¬ 
tion Qi > 2 is respected. In order to verify whether this 
result is particular to this network topology we have an¬ 
alyzed six other multicyclic networks, which do not share 
any particular symmetry [s^. For all these networks, our 
uncertainty relation Qa > 2 is fulfilled. In all cases nu¬ 
merical minimization of Qa leads to a minimum compat¬ 


ible with 2 that is reached in the linear response regime. 
Based on these results, we conjecture that Qq > 2 also 
for general multicyclic networks beyond linear response. 

In conclusion, for nonequilibrium stationary states we 
have conjectured the fundamental limit 2kBT /e^ on the 
minimal dissipation required to generate an output with 
small relative uncertainty e. This bound can be saturated 
close to equilibrium with only one independent affinity 
driving the process. This uncertainty relation provides 
a universal link between the “precision” of a molecular 
machine and the cost of maintaining it. As one specific 
application we have shown how a bound on the num¬ 
ber of states involved in the enzymatic cycle related to 
the Fano factor can be improved provided the affinity 
is known. More broadly, one can expect similar signa¬ 
tures of our fundamental relation quantifying the min¬ 
imal energetic cost for reaching small uncertainty, i.e., 
high precision, to show up in any biochemical or biophys¬ 
ical process at fixed temperature that can be described by 
a Markov network [4l| . Exploring whether and how this 
balance between fluctuations and energetics has guided 
the evolution of chemical reaction networks in living sys¬ 
tems constitutes one intriguing perspective of our ap¬ 
proach. Our fundamental relationship between minimal 
dissipated heat and uncertainty is based, first, on exact 
results in the linear response regime. Second, we have 
proved it for unicyclic networks arbitrarily far from equi¬ 
librium. Third, for multicyclic networks far from equi¬ 
librium we have numerical evidence for several different 
networks. We could not provide a formal proof for arbi¬ 
trary networks, and we expect that the method used for 
unicyclic networks cannot be generalized to multicyclic 
networks, as it requires an expression for the diffusion 
coefficient in terms of the transition rates. 

On the technical level, investigating possible general¬ 
izations of the affinity dependent bound on the Fano fac¬ 
tor in Eq. (HU) to multicyclic networks could lead to 
further new bounds in statistical kinetics. Likewise, it 
would be interesting to explore whether one can derive 
bounds involving higher order cumulants. Finally, we 
emphasize that an algebraic proof of the uncertainty re¬ 
lation in the multicyclic case beyond the linear response 
regime looks like a serious challenge. 
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I. PROVING THAT THE MATRIX G IS POSITIVE SEMI-DEFINITE 


We prove that the matrix G in Eq. (13) in the main text is positive semi-definite. For a given n by n Onsager 
matrix L, where n is the number of independent affinities, we define the (n — l)-dimensional matrix G through 


G — LiiLp-y L\^L\^ 


(25) 


where /3 ,7 = 2,..., n. For convenience and without loss of generality we have set a = 1 in comparison with the main 
text. Since L is positive semi-definite, G is symmetric and its diagonal elements are positive, which are necessary 
conditions for a positive definite matrix. 

From Sylvester’s identity for bordered determinants it follows that 

det(G) =L”f^det(L). (26) 


This relation guarantees that det(G) > 0. A principal minor of G is given by the determinant of the matrix that is 
obtained by eliminating columns and lines from G. For example, a principal minor of order I, denoted by is given 
by the determinant of the matrix that is obtained by deleting column and line /3 from matrix G, with /3 = 2,..., n. 
If we denote by a principal minor of order one of the matrix L, then, similarly to Eq. (ESI), Sylvester’s identity 
for bordered determinants gives 


„( 1 ) _ Tn—3jG) 

9p — ^11 ‘/3 • 


(27) 


Hence, we have that the principal minors of order one of the matrix G are non-negative. Clearly, relations similar to 
(HZl) are valid for any principal minor of order k, of which there are (n — l)!/[fc!(n — 1 — fc)!]. Therefore, all principal 
minors of G are non-negative which, from Sylvester’s criterion, implies G being positive semi-definite. 


II. EXPRESSIONS FOR VELOCITY AND DIFFUSION CONSTANT 


In this section we derive the expressions for velocity and diffusion obtained by Koza 
G“(z) is a iV by TV matrix with elements 


|. The modified generator 


r fcy exp(zd“) 

I - Ej if * = j 


(28) 


It can be shown that the maximum eig envalue of this generator X{z) gives the scale cumulant generating function 
related to the random variable Xa [37|. The current and diffusion constant can be obtained from A( 0 ) with the 
formulas 


J„ = A' (29) 

and 

Dc. = A"/2, (30) 

where the primes denote derivatives taken at z = 0. Furthermore, the characteristic polynomial associated with C°‘{z) 

reads 

N 

P{z,y) =det( 2 //-£“(z)) = ^C„{z)y^. (31) 

n—0 

Since X{z) is a root of this characteristic polynomial, the following relation holds 

N 

^G„(z)A"(z)=0. 

n=0 


(32) 
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Moreover, for z = 0 the matrix in Eq. (l28l) becomes an stochastic matrix leading to A(0) = 0. Taking the derivative 
with respect to z and setting z = 0 in Eq. (15^ leads to 

A' = -C'/Ci. (33) 

Taking a second derivative of Eq. (15^ with respect to z and setting z = 0 gives 

A" = -[C^ + 2C{X' + 2C2{X'f]/{Ci). (34) 

Using Eqs. (1^ . (1501) . (1551) . and (1551) . we obtain the final formulas 

Jc = -C'o/Ci, (35) 


and 


Do. = -{C'o + 2C[J^ + 2C2 J^)/(2C7i), 


which are Eqs. (16) and (17) in the main part. 


(36) 


III. CALCULATIONS FOR UNICYCLIC NETWORK 


For a unicyclic network, the coefficient < 70 ( 2 :) in Eq. (|3T|) is given by the determinant of the matrix (I28L which 
reads Co{z) = —(e^ — l)r_(e-^ — e“*). Therefore, with Eqs. (1551) . (1501) . and the entropy production for the unicyclic 
model a = JA, the quantity Q becomes 


Q = 


2DA 

J 


e-^ + 1 
e-^ — 1 


A-2 


CaT. 


(e-^ 


1 )A 


(37) 


where we used the fact that C( = 0 for a unicyclic network [s^. We now show that C 2 T-/C 1 reaches its maximum 
when = (r+)^/^ and k~ = (T”)^/^ for i = 0,..., — 1. It is convenient to write the rates in the form 


kt = 




and 


k, = 


,-Aej2 


(38) 


where the constraint 


1 fixes the affinity A. We define the function 


where 


h{n,s) = 


0(n, s) 


ELn+i - Er=o if s > n + 1 

0 . + E:=o - EL .+1 0 . if s < n. 


(39) 


(40) 


and 


N-l 




(41) 


2=0 


With this function the coefficients of the characteristic polynomial in Eq. (ED) can be written as 

AT-l N-l 


c”! = E E 


n—O s—0 


and 


N-lN-l-l n+Z 

E E E 

1 — 1 n—0 s=n+l 


N-l 


C 2 T-=Y^ E E Hn,s)('^h{n + l,t)+ Hn + l,t) 


i=0 


t—n-\-l-\-l 


(42) 


(43) 





These two relations can be obtained from the formula 


N-l 

Ci=Y, Y[ {k+r-{k-f-Ma,b) (44) 

a,b m,n=0 

where Z = 1, 2, a (6) denotes a vector with components a„ € {0,1} {hn € {0,1}) and tpi{a, b) € {0,1}. The function 
ipi{a, b) is non-zero only if + bn) = N — I and if for all n = 0,1,..., TV — 1 the relations a„ = &„ = 1 and 

«n = bn +1 = 1 are not fulfilled. Moreover, a graphical representation of Eq. (1421) can be found in [s^ . 

To find the maximum of f{{di}, {4>i}) = C- 2 V-/C\ we consider the function 




N-l 


1-E 


(45) 


\ i=0 / 

where Aq is the Lagrange multiplier. This function is maximized for {&*}, {4’i}^ which are given by the solution of 
the equations 


d^jm},m) = o (46) 

and 

dej{m,m}) = ^0, (47) 

for j = 0,1,..., — 1. Before taking derivatives of the function / it is convenient to define 

hj{nj,Sj) = h(n,s), (48) 

where rij = {n — j + N) mod N, sj = {s — j + N) mod N, and j = 0, 1, ..., — 1. Due to the symmetry of the 

unicyclic network, we can rewrite Eqs (1421) and (1431) in the forms 

7V-1 N-l 

Cl = ^ ^ hj (n,s), (49) 

n—0 s=0 


and 


N—IN—1 — 1 n+Z / n N-l 

CsT. = ^ E E bj{n,s)('^hj{n + l,t)+ ^ hj{n + l,t) 

I—I n—0 s—n+1 ^ t—0 

which are valid for all j. A first derivative of the function / reads 

^ ^ ^ a.(C2r_) _ 2(a,Ci)C2r_ 


c? 


c? 


where x = (j^j^Oj. Taking derivatives of Eqs. (|49l) and (l5Ql) with respect to (j)j, leads to 


N-l 


= -[ClhMo) 


n—0 


and 


a..C2r_ = ^ 


N-lN-l-1 n+Z 

2 C 2 r_-^ ^ ^ h,{n,s)h,{n + l,0) 

1—1 n—0 s=n+l 


(50) 


(51) 


(52) 


(53) 


At the symmetric point 6* = 1/N and (j)* = (j) for all i, these derivatives become 
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and 

d^,C2T. = ^^C2r- (55) 

Using these two last relations in Eq. m we obtain that d^.f — 0 at 0* = 1/N and cj}* = 4>. Taking derivatives of 
Eqs. ((4^ and (l50l) with respect to 9j leads to 

. Af-l / n N-1 

dg^ C*! = y X! ( XI “ X 

n—0 \s=l s=n+l 

and 

n+/ / N-1 

dg^C 2 T- = --Y^ Y X hj{n,s)i2 Y hj{n + l,t) + hj{n +1,0)^ (57) 

1—1 n—0 s—n+1 \ t—n-\-l-\-l / 



From these two relations and Eq. (ICT) we obtain that the derivative dg.f is independent of j for 0* = 1/fV and 4>* = (p 
for all i. In this case the Lagrange multiplier in Eq. (iTTll becomes 


An — 


AA 


2(e-A/W-l)2(e-A_i)3]v2’ 


(58) 


where 


A = 2e-^/" - e2(^+i)-^/^(7v - 1)N - e‘^^N{N + 1) - N{3N + 1) + 

e-^(A^ - SiV^) + ei^N+i)A/N^{ n+i)a/n 

Hence, we have shown that for symmetric rates, with 9* = 1/N and </>* = </) for all i, Eqs. (l46l) and (iTfl) are satisfied. 
From the explicit calculations of the derivatives they seem to be the unique solution, however, we are not able to 
provide a rigorous proof of this uniqueness. The maximal value of the function / is 

C 2 T- (e^+^/^-l)(7V-l)-(e-^-e-^/^)(7V + l) 

~Cf~ ~ 2N{e^/N - l)(e^ -1)^ ' 

Inserting this expression in Eq. (1571) leads to the lower bound in Eq. (20) in the main text. 


IV. VERIFYING THE BOUND Q > 2 FOR MULTICYCLIC NETWORKS 


We present numerical results supporting the uncertainty relation for six different models. The simplest example of 
a multicyclic network is shown in Fig. 2(a) The two states are denoted by I and 2 and the transition rates between 


then are k '/2 and fc^i, with v = a, b,c for the three different links. We take the output to be the number of transitions 
from 1 to 2 minus the number of transitions from 2 to I in link a. With this choice the modified generator in Eq. 
becomes 


— (h°- 4- _l_ t-c t l-a pZ , Lb i Lc 

V"'12 ' "'12 ' "'121 "12'= W "12 w "12 

La p-z I 1,6 I jUc 4- lA 4- \ 

"21'= w "21 w "21 l"21 w "21 w "21/ 


(61) 


The related current Jq and diffusion coefficient Da are obtained from Eqs. 


and (1551) . They read 


and 


for this model reads 


■ K2)/Ks 

(62) 

- (Kl - K2f/{KN: 

(63) 

1 Ic^ —1— 1 Jc^ 1 

^12 ^ ^12 ^ ^12 ^ ^21 ^ ^21 ^ ^21- 

The entropy production 


cr = Ja In 


La Lb 
'^12'^2i 
ua ub 

ft-oi Avi r 


Jc In 


^12^2 
Avoi Ail o 


(64) 
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FIG. 2. Network of states for multicyclic models, (a) Two-state model with two fundamental cycles, (b) six-state model with 
two fundamental cycles, (c) five-state model with three fundamental cycles, (d) four-state model with two fundamental cycles, 
(e) fully-connected four-state model, (f) three-state model with four fundamental cycles. 


where J^, = Pik '(2 — P 2 ky^ ■ W ith Eqs. (15^ . (1551) . and (IMl) . we obtain Qa = 2DaCr/Ja- A numerical test of the bound 
Qa > 2 is given in Fig. 3(a)[ where we evaluate Qa for randomly chosen transition rates. 

Following the same procedure explained above for the two-state model we calculated Q for the other five network 
of states in Fig. [2] Numerical tests of the inequality Qa > 2 for these networks are shown in Figs. |31 with the 
chosen observable Xa for each network indicated in the caption. Besides evaluating Qa at randomly chosen transition 
rates we have also minimized it numerically. In all cases the minimum value is compatible with 2. Furthermore, the 
minimum is reached at the linear response regime with the affinity Aa small and the other affinities even smaller, 
com patib le with A/? = 0 for /3 a. We point out that for the symmetric networks in Fig. 2(a)[ Fig. |2(e)[ and 
Fig. 2(f) I all the currents are equivalent, i.e., all links connect the same number of nodes and all nodes have the same 
number of links. For the other networks, we have also checked the bound for the following cases: outputs X 13 and 
X34 for the networ k in F ig. |2(b)| outputs X41, X12, and X34 for the network in Fig. 2(c)[ outputs X^^ and A 23 for 
the network in Fig. |2(d)[ 
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log i,, 

(d) 






log I<^'i2 
(f) 


FIG. 3. Numerical check of the bound Q> 2. (a) Network of Fig. 2(a) with output X 12 , (b) network of Fig. 2(b) 
X 12 , (c) network of Fig. 2(c) with output X 23 , (d) network of Fig. 2(d) with output X 12 , (e) network of Fig. 2(e) 


with output 
with output 

X 12 , and (f) network of Fig. 2(f) with output X'^ 2 - In nil cases the transition rates kij are randomly chosen by generating a 
random number x between —3 and 3 and then taking kij = 10^. 




































